Reversible Transport of Interacting Brownian Ratchets 
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The transport of interacting Brownian particles in a periodic asymmetric (ratchet) substrate is 
studied numerically. In a zero-temperature regime, the system behaves as a reversible step motor, 
undergoing multiple sign reversals of the particle current as any of the following parameters are 
varied: the pinning potential parameters, the particle occupation number, and the excitation am- 
plitude. The reversals are induced by successive changes in the symmetry of the effective ratchet 
potential produced by the substrate and the fraction of particles which are effectively pinned. At 
high temperatures and low frequencies, thermal noise assists delocalization of the pinned particles, 
rendering the system to recover net motion along the gentler direction of the substrate potential. 
The joint effect of high temperature and high frequency, on the other hand, induces an additional 
current inversion, this time favoring motion along the direction where the ratchet potential is steeper. 
The dependence of these properties on the ratchet parameters and particle density is analyzed in 
detail. 

PACS numbers: 05.40.-a, 87.16.-b 



I. INTRODUCTION 

Increasing interest in the manipulation of small objects 
has attracted much attention to the so called Brownian 
ratchets [1 0, i, i, Hi- These systems are able to pro- 
mote the directed transport of particles on a periodic 
substrate lacking reflection symmetry by means of sym- 
metrical, unbiased, non-equilibrium fluctuations. Their 
applicability can range from understanding the molecu- 
lar transport of proteins 0] to the controlled transport 
of flux quanta in superconductors [1, [13, [ill. IT^ and 
separation of mixtures in nanometric devices [l3[ . Albeit 
equilibrium (thermal) fluctuations alone cannot induce 
such a motion, it plays a key role within transport effi- 
ciency when the system is driven out of equilibrium. 

Although these systems are designed with a nominal 
"easy" direction of motion (a gentler slope than that of 
the ratchet potential), in many cases it is not a trivial 
task to predict even the real drift direction. In a seminal 
paper, Doering illustrated that when the system is 
driven by a fast non-equilibrium (but stochastic) fluctu- 
ating force, the particle current may be reversed towards 
the "hard" direction, depending upon the noise statis- 
tics. Much later, a similar effect was observed for ratch- 
ets driven by a high frequency harmonic force [^ . In this 
later case, the reversed motion resulted from the inter- 
play between fast driving and equilibrium white noise, as 
demonstrated analytically by Fistul [l3l- In both situa- 
tions described above, the current reversal was induced 
by stochastic processes. 

When several particles arc put in a ratchet substrate, 
their mutual interaction c an g ive rise to a variety of co- 
operative phenomena [TBI . Il6l | Remarkably, even when 
the substrate potential is perfectly symmetric, directed 
macroscopic transport of the interacting particles is pos- 
sible. This can be accomplished, for instance, by spon- 
taneous symmetry breaking when the system is driven 
out of equilibrium [T^ or by inducing correlated fluctua- 



tions in the interparticle bonds |17| . Interparticle inter- 



actions can also induce deterministic current reversals. 
A dramatic example is the recently discovered multiple 
drift reversals of interacting ratchets [l2|- Such an ef- 
fect was observed experimentally in a system of vortices 
captured by a nanocngincered periodic array of asymmet- 
ric pinning sites in a superconducting film. The vortex 
drift direction was seen to change sign several times as 
a function of the vortex occupation number. Numerical 
simulations suggested that such a sequence of reversals 
is quite general and results from the interplay between 
particle-particle and particle-substrate interactions. In- 
deed, other systems of interacting ratchets seem to ex- 
hibit similar behavior [l^, [l^, [13] ■ This reversed current 
phenomenon has brought about a high level of control- 
lability over the motion of small particles. Although the 
necessary ingredients seem to be clearly set, the dynami- 
cal mechanism behind this effect is not completely under- 
stood. Furthermore, the work carried out so far has not 
explored the influence of thermal noise and non-adiabatic 
oscillating forces. 

In this paper, wc set out to illustrate how the drift of 
long-range interacting small particles on a ratchet sub- 
strate can be manipulated by suitably controlling the in- 
terplay between particle-particle and particle-substrate 
interactions, thermal noise, excitation amplitude, and 
frequency. As demonstrated in Ref . [l2| , the drift of re- 
pelling particles undergoes multiple sign reversals, from 
positive (along the easy ratchet direction) to negative 
(along the hard ratchet direction) as the number of par- 
ticles per ratchet period (n) changes from odd to even 
for low amplitude of the external oscillating force and no 
thermal noise. In the present work, we show that, at a 
fixed particle density (n > 2), a sequence of n — 1 re- 
versals emerges as the excitation amplitude is increased 
above the ratchet threshold. By detailed characterization 
of the dynamical states, wc conclude that such rever- 
sals are induced by dynamical transitions of the particle 
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FIG. 1: Double-sine potential used in the simulations (Eq.[2|, 
for /? = 0.25, 0.5, 1.0 and 1.5. The curves are vertically shifted 
for clarity. Energy in units of the interparticle interaction 
energy scale, Eq, and position in units of the ratchet period. 



chain. Furthermore, wc undertake a detailed analysis of 
how thermal fluctuations assist the restoration of the pos- 
itive current and how the combination of thermal noise 
and a high excitation frequency leads ultimately to the 
negative ratchet drift. 



II. SYSTEM AND SIMULATIONS 

We consider interacting particles confined in a one- 
dimensional asymmetric periodic substrate potential and 
subjected to a fluctuating driving force. We assumed that 
the particles evolve in time according to the Langevin 
equation of motion: 

mv, = -9,. ([/, + U^J) - IV, + Fd{t) + (1) 



where 7 represents the viscosity coefficient, m is the par- 
ticle mass. Us is the substrate potential, Uij is the pair 
potential between particles i and j, Fd{t) is the fluctuat- 
ing driving force and (,i{t) is the uncorrelated Gaussian 
noise, which accounts for the thermal fluctuations at a 
substrate temperature T via the fluctuation-dissipation 
theorem: {£.i{t)ij(t')) = 2-fkBTS,jS{t - t'). 

For the substrate potential, C/s^we assumed the usual 
asymmetric double-sine function [5|, 



Here Uso and (3 arc free parameters controlling the in- 
tensity, the asymmetry and the internal structure of the 
ratchet potential. Some plots of Uso{x) are shown in 
Fig. [1] The shape of the potentials can range from one 
with a single local minimum, for (3 < 0.5, to one with two 
local minima per period, for (3 > 0.5. The interaction po- 
tential between particles i and j is given by a repulsive 



Such an interaction is relevant for example, for vortices 
in superconducting thin films. The external driving force 
is given by Fd = Asinicut), where A and tu are the ampli- 
tude and the frequency respectively. Hereafter, we shall 
adopt the following unit system: Eq, for energy, a, for 
length, to = ^/ma^/Eo, for time, and ujq ~ Stt/^o, for 
angular frequency. 

The simulations were performed on a section of the 
substrate with size i = 10 and a total number of particles 
N ranging from 10 to 40 (n = 1,2,3 and 4) under peri- 
odic boundary conditions. The simulations were carried 
out using a standard finite difference method adapted for 
stochastic equations. For a given set of parameters, we 
discarded the first 10^ to 10^ time steps to allow the sys- 
tem to relax to the dynamical steady state and follow 
up the same number of time steps for the time averages. 
To enhance the signal-to-noise ratio, averages were cal- 
culated of the physical quantities for several realizations 
of the thermal noise (up to 50 samples, for low excitation 
frequencies, and up to 700 samples for high frequencies). 
The mean particle current, (u), is calculated by taking 
the average of the particle center-of-mass (cm.) veloc- 
ities, Wcm = over time and over the thermal 
noise realizations. To characterize the dynamical states, 
the time-averaged density function p(x) ~ S{xi)) and 
the cumulative mass m{x) = p[x')dx' were calculated. 

To set a few parameters, we chose m = 1 and 7 = 16 
in most simulations. Other values of m and a compari- 
son between underdampcd and ovcrdamped dynamics is 
discussed in Sec. IHI Dl 



III. RESULTS 

The deterministic (zero temperature) phase diagrams 
shown in Fig. [2] indicate the preferred direction of the 
mean particle current for a given occupation number n, 
ranging from 1 to 4, as a function of Uso and /3, in the 
limit where the applied force triggers directed motion, 
i.e. in the ratchet threshold. This measurement was 
performed by adiabatically increasing an external force 
until a collective motion of the particles was observed. 
The measured threshold forces for motion to the right, 
Fr, and to the left, Fi, were then used to compute the 
effective asymmetry parameter a = 1 — Fr/Fi, which is 
positive (negative) for rightward (leftward) motion. The 
diagrams exhibit the usual sign reversals as a function of 
n and of the ratchet potential parameters Usq and (3 (l^ . 



A. The zero-temperature adiabatic regime 

Here, the ratchet dynamics is explored when the driv- 
ing force amplitude A is increased within the zero- 
temperature, low-frequency (adiabatic) regime. Fig. [3] 
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FIG. 2: (Color online) Density plots of the (dimensionless) 
effective asymmetry a = 1 — Fr/Fi as a function of the ratchet 
parameters UsO (in units of -Eo) and j3 (dimensionless), for n = 
1 to 4. Fr (Fi) is the critical adiabatic depinning force applied 
to right (left) in the sample. For f3 > 0.5 and Uso/Eq < 1.0, 
the current direction is positive for odd n and negative for 
even n. 



shows the particle current function of ^ for n = 1, 
2, 3, and 4. The calculations were performed for the 
ratchet potential parameters Ugo ~ 0.74 and f3 = 1.5 
(upper curve in Fig.[T]) and a frequency of u; = 0.01. The 
deterministic (T = 0) solutions are given by the solid 
lines. The noisy adiabatic regime will be discussed fur- 
ther on. For n > 1, (n — 1) sign reversals occur with 
A increasing. Such reversals take place well above the 
ratchet threshold, contrary to those observed in Fig. ^ 
and therefore are induced by dynamical transitions in the 
array of particles. 

In order to conduct a detailed analysis of the possi- 
ble dynamical states and dynamical transitions and how 
they lead to the observed sign reversals, the dynamics 
of particles driven by constant force values was simu- 
lated. The force intensity was slowly ramped up with 
the force applied to the positive and then to the nega- 
tive direction. In this case, the net current is given by 
J{F) — J+{F) + J_(F), with J+(-) being the mean par- 
ticle current for the force applied to the right (left) . Such 
dc measurements of the ratchet effect is equivalent to an 
ac measurement where the oscillating force is an adia- 
batic square wave [l^. The advantage of such a strat- 
egy is that the dynamical states can be characterized for 
a constant force during an arbitrary long time interval, 
thus allowing one to compute statistical quantities for 
that specific force value. 

The net current J, calculated for the dc measurement 
described above, as a function of force intensity F for 
UsO = 0.74, (3 = 1.5, and n = 3, is presented in Fig. H) 
The J{F) curve qualitatively resembles the J{A) curve 
of Fig.[5]-(c), calculated for sinusoidal drives. 

Fig. m presents plots of the particle density p{x) (panel 
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FIG. 3: Net particle current J as a function of amplitude A of 
the sinusoidal drive for the ratchet potential parameters Uso = 
0.74 and /3 = 1.5. Full lines indicate T = and line-symbols 
indicate T = 0.625, in units of Eo/ks- For n = 1 (a) there is 
no current inversion whereas for n — 2 (b), n — 3 (c) and n — 
4 (b) respectively, there are clearly n — 1 current inversions 
induced by amplitude for the zero-temperature case. 
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FIG. 4: (Color online) Net particle current J as a function of 
force intensity calculated for a dc measurement of the ratchet 
effect (see text). Arrows indicate the force intensities used in 
Fig.[S] 



(a)) and of the cumulative mass m(x) (panel (b)) for 
n = 3 and positive driving forces with intensities F = 9, 
17.5, and 22. As indicated in Fig. IH these force intensi- 
ties correspond to low-amplitude positive current, nega- 
tive current, and high- amplitude positive current, respec- 
tively. The region of low-amplitude positive drift is char- 
acterized by a double peak at each local minimum of the 
substrate potential, which is represented in Fig. [5] as the 
background area plots. The jumps in m(a;), integrated 
for F = 9, is evidence that each double peak observed in 
p[x) corresponds to exactly one particle distributed in the 
close vicinity of the substrate local minima, i.e. two lo- 
calized particles per substrate period, while the smoother 
regions of mix) correspond to a third particle spread out 
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across each substrate period, i.e. one delocalized particle 
per substrate period. This leads to the conclusion that 
all minima of the substrate potential arc occupied most 
of the time, except for a very short time interval during 
which a particle trapped in one minimum is knocked out 
and replace by another. Therefore, the effective dynamics 
may be roughly described as N/3 (delocalized) particles 
moving over an effective ratchet potential composed of 
the bare substrate potential Us plus the potential due 
to the other 2N/3 (localized) particles fixed at the local 
minima of the substrate potential. Clearly, this effective 
potential preserves the symmetry of the substrate po- 
tential, favoring the observed positive current sign. For 
F = 17.5, however, the p{x) peaks at the weaker minima 
are drastically reduced and the mass distributes broadly 
near these minima in comparison to the sharp distribu- 
tion near strong minima. Now. the effective potential 
comprises particles localized at the stronger min- 
ima, rendering an inversion of the effective asymmetry of 
the system. Finally, at = 22, all p{x) peaks are small 
and there is no point in the substrate potential where 
m{x) is sharply distributed. Hence all particles are now 
delocalized and the system restores its original symmetry. 

For the negative drive direction, a similar sequence of 
dynamical transitions was observed, but with different 
critical values of F. An overall picture of what happens at 
positive and negative drive directions and the relevance 
for the observed ratchet states is given in Fig. [5] where 
the main peaks located at the deeper (p(.Ti)) and the 
shallower {p{x2)) local minima arc plotted against force 
intensity for different drive orientations. Roughly speak- 
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FIG. 5: (Color online) Time-averaged particle density p[x) 
(a) and (dimensionless) cumulative mass m{x) (b) for posi- 
tive, constant drive values F = 9 (full lines), 17.5 (dashes), 
and 22 (dash-dots). The background area plots represent the 
substrate potential used in this calculation. 
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FIG. 6: (Color online) Local maximum values of the time- 
averaged particle density as a function of force intensity for 
positive (top) and negative (bottom) drive directions. The 
maxima were measured in the vicinity of the strong (thick 
lines) and the weak (thin line) potential wells (see Fig. [5)l . 



ing, the negative current region is characterized most of 
the time by only one localized particle per substrate pe- 
riod, which tends to invert the system symmetry inducing 
the negative current, while the positive current regions 
are mostly characterized by either two (low amplitude) 
or no (high amplitude) localized particles, which tends 
to keep the system symmetry. The exact force intensi- 
ties where the sign inversions take place arc a non-trivial 
result from the competition between the dynamics at pos- 
itive and negative drive direction. 

More generally, at small drive amplitudes, for any 
n > 1, there are (n — 1) particles in a substrate period 
that can be considered to be localized whereas one par- 
ticle is delocalized. As the amplitude increases, the lo- 
calized particles are dynamically delocalized one by one, 
inducing a total of {n — 1) sign reversals. In other words, 
the system is sequentially reduced to dynamical states 
similar to those observed at the small amplitude regime 
of lower particle density systems. 



B. Noise-induced sign inversions 

To check the effect of thermal fluctuations in the adi- 
abatic limit, detailed simulations of Eq. [1] have been 
performed for T > and low frequency [uj — 0.01). 
Fig. [3] presents the mean particle current calculated for 
T = 0.625 (line-dot), which corresponds to the strong 
noise limit. For this temperature value, the ratchet cur- 
rent is essentially positive for all n, particularly in the 
amplitude ranges where the current direction at T = 
is negative, that is, for such amplitude ranges, thermal 
noise induces a current reversal. 

To investigate these noise- induced inversions in greater 
detail, we calculated the temperature dependence of the 
net current, J{T), of the n — 2 particle array for 
UJ = 0.01, P = 1.0, and several values of the substrate 
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strength Ugo and the drive amphtude A. Fig. [7] (a) maps 
the region in the Usq-A plane where the T = drift state 
is negative. The rectangle depicts the region considered 
for the T > simulations. The results of the J(T) sweeps 
were used to build the T-Usq-A diagram shown in Fig. [7] 
(b) , where the surface represents the transition tempera- 
ture, Ti„y(A, Uso), between the negative and the positive 
drift states. 

Interestingly, T!i„„(A) seems to decrease linearly with 
the amplitude. In other words, the critical amplitude, 
Ac, at which the thermally- induced drift inversion oc- 
curs, scales as Ac = a — bT. This becomes evident in 
Fig.[7]-(c) where Ac{T) is plotted for several values of UsO. 
To understand this behavior, we first recall the results of 
the previous paragraphs, which show that the inverted 
ratchet drift observed for a n = 2 particle array can be 
described as half of the particles moving over the poten- 
tial generated by the substrate and the remaining (local- 
ized) particles. Therefore, one can roughly associate Ac 
with the force where such dynamics is destroyed, that 
is, where the localized particles delocalize, thus restor- 
ing the original substrate symmetry. On the other hand, 
the linear behavior of Ac{T) suggests a transition to a 
regime of creep motion resulting from thermal activa- 
tion of the localized particles over the energy barriers of 
the (tilted) substrate potential. In systems exhibiting 
similar dynamics, such as flux lines in hard supercon- 
ductors [2l| and domain walls in ferromagnetic materi- 
als [l^l, the critical force for a detectable creep motion 
decreases linearly with T as long as the relevant energy 
barrier at zero drive, AUs, has no explicit temperature 
dependence and fcsT ^ AUg H^. In the present case, 
for (3 = 1.0, AUs ~ 4[/so, which is much bigger than ksT 
near the transition surface, as can easily be inferred from 
Fig.EKb). 



placement of a/n per drive cycle. This result provides 
the opportunity for controlling the particle current ac- 
curately, at a given frequency, by changing the particle 
density. 

At non-zero low temperatures the lock-in effect is 
smeared out by thermal fluctuations. At sufficiently 
high temperatures, thermally induced current rever- 
sals, similar to those observed in the adiabatic regime, 
are observed in the amplitude range where the zero- 
temperature drift direction is negative (Fig.[8ja)). How- 
ever, when the frequency is even higher, the combina- 
tion of fast rocking (lu ~ 7) and thermal fluctuations 
(T = 0.625) leads to another reversal, resulting in the 
negative current sign in the full amplitude range. This re- 
sult is somewhat expected since, as mentioned above, at 
a high T the particles may be considered independent of 
one another and, for non-interacting rocking ratchets at 
high temperatures, a sufficientl y h igh driving frequency 
induces the negative current fsl. Il4l|. 

Fig. [9] illustrates the sequence of sign inversions when 
one first increases the temperature in the adiabatic (uj = 
0.01) regime (a) and then, at a fixed high temperature 
value (T = 0.625), increases the frequency up to w = 10 
(b). Note that the current induced by high frequencies 
is approximately one order of magnitude smaller than 
that induced by temperature. Indeed, the frequency- 
induced reversed current of non-interacting ratchets is 
usually small because such inversions occur at frequencies 
much higher than the natural frequency of the system. 
The same behavior is observed here because at a high 
frequency and temperature the system is essentially non- 
interacting. Interacting ratchets in the adiabatic regime, 
on the other hand, offer several ways to reverse the sign 
of the current (by changing either the particle density, 
the force amplitude or the temperature) while keeping 
its absolute value at the same order of magnitude. 



C. The high frequency regime 

Here we turn our attention to the effects of high fre- 
quency drives. To fix some quantities, we set n — 2, 
UsO = 0.74, and = 1.0. The amplitude dependence of 
the mean current for a few frequency and temperature 
values is shown in Fig. [51 

At T = 0, J{A) is characterized by a series of plateaus 
and spikes produced by the well-known lock-in effect ob- 
served in single-particle rocking ratchets. Notice, how- 
ever, that the mean particle velocity assumes discretized 
values at multiples of vq = 0.5wa/27r = 0.016, which 
means that each particle moves across an average dis- 
tance of a/2 during a cycle of the external drive. This 
result is in contrast with that of single-particle ratch- 
ets, where the lock-in state corresponds to the particle 
moving a step forward across a full period of the sub- 
strate in one drive cycle, which gives wq = loajl'E Q. By 
performing similar simulations for other particle occupa- 
tion number values, we have observed that, in general, 
vq = Loajl'Kn, corresponding to an average particle dis- 



D. Overdamped versus underdamped dynamics 

In the previous paragraphs, we explored the dynamical 
properties of ID arrays of particles of unity mass moving 
on a ratchet substrate with friction coefficient 7 = 16. 
Here, we shall discuss other values of m, spanning the 
overdamped and underdamped regimes. Although the 
ratchet threshold diagrams calculated for the different m 
values are very similar to those presented in Fig. [5] for 
TO = 1, the behavior of the moving phases are qualita- 
tively different for different intervals of to values. 

As noticed by Borromeo et al. [l^l , the particle current 
(for rocking, non-interacting ratchets) behaves qualita- 
tively different from the usual overdamped (to = 0) limit 
when the ratio 7/771 becomes smaller than the libration 
frequency of th e par ticles at the bottom of the ratchet 
wells (wiib = yjk/m, with k the spring constant of the 
well). In this moderate underdamped regime, the onset 
of the ratchet current becomes abrupt and the ratchet 
tail in {J) {A) decays faster as compared to the {J){A) 
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FIG. 7: (Color online)(a) Zero-temperature phase diagram of the drift direction of a n = 2 chain for /3 = 1.0, — 0.01. (b) 
3D phase diagram of the same system showing the inversion temperature surface separating the negative (low temperature) 
and positive (high temperature) current regions eis a function of Uso and A. The ranges of Uao and A are represented by the 
rectangle in (a), (c) Temperature dependence of the critical amplitude, Ac, where thermally assisted inversions take place, 
calculated for several Uso- 



of an overdamped ratchet. When 7/m ^ i^iib, however, 
the above described behavior is substituted by a scenario 
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FIG. 8; Mean particle current of the n — 2 chain as a function 
of amplitude for (a) moderate frequency, u — 0.2, and (b) 
high frequency, u = 7.0. Substrate parameters are (3 — 1.0 
and Uso = 0.74. Full lines indicate T = and symbol-line 
plots indicate T = 0.5 (a) and T = 0.625 (b) in units of 
Eo/kB. 



dominated by chaotic dynamics. 

In our case, with the ratchet potential given by Eq. [2] 
and choosing /3 — 1.5, one obtains k = 262Uso/a^, so 
that the critical mass that roughly separates the under- 
damped and overdamped regimes is m* = j'^/k = 1.32. 
Fig. [TOl shows the (J) (A) characteristics for particle den- 
sities 71 = 2, 3 and 4, and different mass values, ranging 
from to 5. Larger m values lead to chaotic dynam- 
ics, which is out of the scope the present paper (in fact, 
already for m = 5.0, chaotic trajectories were observed 
for a few drive amplitude values resulting on the fluc- 
tuations exhibited by the (J) (A) curves for m = 5.0 
and n = 4). As it is clear, for m > m* (i.e. m — 2.0 
and m = 5.0), the transitions between different ratchet 
phases are more abrupt than those observed for m < 1.0. 
Notwithstanding, the n — 1 oscillations are still present 
and the behavior of {J) {A) for m = 2.0 and m = 5.0 is 
quite revealing with respect to the scenario described in 
Sec. nil Al For instance, the interval A > 8.7 in the 
n — 4, m ~ 5.0 curve seems to qualitatively imitate 
the whole n = 3, m = 5.0 curve, whereas the interval 
A > 12.3 in the n = 3, m = 5.0 curve seems to qualita- 
tively imitate the {n ~ 2,m = 5.0) curve. This is consis- 
tent with the fact that it is the number of localized par- 
ticles which determines the ratchet direction. The qual- 
itative similarity between curves with different n arises 
from the fact that in the underdampcd regime the ratio 
between the particle current and the force amplitude has 
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FIG. 9: (Color online) (a) Mean particle current, J, of the 
n = 2 chain as a function of temperature, T, in the adiabatic 
regime {lo = 0.01) for /3 = 1.0, UsO = 0.74, and A = 10. (b) 
Frequency dependence of J for the same substrate parameters 
and A = IQ starting from the point indicated by a circle in 
(a) (r = 0.625). 



a universal shape which does not depend on the details 
of the potential but rather on its critical forces [l^] ■ 

For m < TO*, the curves are much smoother and, for 
n = 3 and n = 4, the high-amplitude negative drift phase 
becomes positive for to < 0.5. Notwithstanding, the os- 
cillations, reminiscent of the interplay between localized 
and dclocalizcd particles, are still present, even in the 
TO = case. 



the symmetry of the effective ratchet potential and, con- 
sequently, the preferential drift direction. In the non- 
zero th ermal noise regime, the oscillating behavior o f 
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FIG. 10: (Color online) Net particle current J as a function 
of amplitude drive A for n = 2 (a), n = 3 (b) and n = 4 
(c) calculated for different mass values: m = (full lines), 
0.33 (dashes), 0.5 (short-dashes), 1.0 (dash-dots), 2.0 (dash- 
double-dots) and 5.0 (long dashes). The ratchet parameters 
are the same as in Fig. [3] 



IV. SUMMARY 

In summary, the transport properties of interacting 
particles confined in a ID ratchet potential has been 
investigated in detail. It has been shown that, in the 
low-temperature adiabatic regime, the particle chain of 
density n undergoes a series of n — 1 current inver- 
sions as the oscillating drive amplitude is increased for 
a range of the investigated parameters. A systematic 
study of the time-averaged particle density and cumula- 
tive mass, calculated at different drive intensities, sug- 
gests that these inversions result from the competition 
between different dynamical states of the chain where a 
certain fraction of the particles are effectively trapped 
by the substrate, while the others float over the effective 
potential generated by the substrate and the localized 
particles. The number of localized particles determines 



the amplitude-dependent mean current observed at low 
temperatures gives place to a positive-drift phase in the 
whole amplitude range at a sufficiently high temperature. 
In cases where the drift is negative at T = 0, the cur- 
rent direction becomes positive at a certain temperature- 
dependent critical amplitude Ac{T), which decreases lin- 
early. Such behavior suggests that these noise-induced 
inversions can be associated with thermal activation of lo- 
calized particles over the substrate barriers. Finally, the 
high-frequency regime has been addressed for both weak 
and strong noises. For weak noise, the steps and spikes 
reminiscent of the well-known lock-in effect have been ob- 
served. But here the size of the steps can be fine-tuned 
by changing the particle density. In the strong noise case, 
the interplay between thermal fluctuations and the fast 
oscillations of the driving force leads to the complete in- 
version towards the negative drift direction in the whole 
amplitude range. 
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